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ABSTRACT 

We delineate the techniques and prospects for using Higgs 
pair production in e + e~ or collisions to probe GUT-scale 

boundary conditions in the minimal supersymmetric standard 
model. 

I. Introduction and Results 

The heavier CP-even, the CP-odd and the charged Higgs 
bosons (H°, A , and H , respectively) of the minimal super- 
symmetric standard model (MSSM) (see Ref. [0] for a recent 
review) are typically predicted to be fairly heavy (with m#o ~ 
m^o ~ m H ± > 200 GeV) in models where electroweak sym- 
metry breaking at scale mz arises as a result of evolution be- 
ginning from simple GUT/Planck-scale boundary conditions. 
Nonetheless, once y/s is large enough that e + e~ — > H°A° and 
e + e~ — » H + H~ (or their analgoues) are kinematically 

possible, event rates are substantial for expected machine lu- 
minosities, and discovery and study of these Higgs bosons be- 
comes possible [0]. The all-jet and high-multiplicity final states 
coming from H", A -> bb, tt and H + -> ib, H~ -> bt are 
background free and for the model we study provide appropri- 
ate and efficient signals with rates that are adequate even when 
SUSY decays are present. Further, in the all-jet channels, the 
individual Higgs boson masses, m^o, m H o and m H +, can be 
measured. Event rates and decay branching fractions are typi- 
cally such that it will be possible to 'tag' one member the pair 
in such a fully reconstructable final state and then study the de- 
cays of the untagged member of the pair. Here, we point out the 
very dramatic sensitivity of measurements of decay branching 
fractions to the GUT boundary condition scenario, illustrating 
in particular the high statistical level at which various not ter- 
ribly different scenarios can be distinguished from one another 
using ratios of branching fractions. A more detailed treatment 
of this analysis appears in Ref. Q. 

In the simplest GUT treatments of the MSSM, soft supersym- 
metry breaking at the GUT scale is specified by three universal 
parameters: 

• mo: the universal soft scalar mass; 

• m 1( / 2 : the universal soft gaugino mass; 

• Aq: the universal soft Yukawa coefficient. 
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The absolute value of /i (the Higgs mixing parameter) is deter- 
mined by requiring that radiative EWSB gives the exact value 
of mz for the experimentally measured value of rat, however, 
the sign of /i remains undetermined. Thus, the remaining pa- 
rameters required to completely fix the model are 

• tan f3: the vacuum expectation value ratio; and 

• sign(/x). 

We remind the reader that a universal gaugino mass at the GUT 
scale implies that M3 : M 2 : Mi ~ 3 : l ; 1/2 at scale ~ mz. 
For models of this class one also finds that [i M± t 2. These 
two facts imply that the xi is mainly bino, while \ 2 and \t are 
mainly wino, with heavier gauginos being mainly higgsino, so 
that m~ n ~ m~+ ~ 2m~ n . 

2 ^1 X-l 

We will consider three representative GUT scenarios char- 
acterized by increasingly large values of too relative to m 1 / 2 
(which translates into increasingly large slepton masses as com- 



pared to to- 



and to~+): 

Xl 



• "No-Scale" (NS): A = m = 0; 

• "Dilaton" (D): m X j 2 = —Aq = y/3mo; 

• "Heavy-Scalar" (HS): too = m 1 / 2 , Aq = 0. 

Within any one of these three scenarios, the model is completely 
specified by values for m 1 / 2 , tan/3andsign(/i). We will present 
results in the (m 1 / 2 ,ta,n (3) parameter space for a given sign(/i) 
and a given choice of scenario. Our notation will be NS~ for 
the no-scale scenario with sign(/i) < 0, and so forth. 
In exploring each of these scenarios, we proceed as follows. 

• First, we delineate the allowed region of (mx/2, tan /3) pa- 
rameter space consistent with all available experimental 
and phenomenological constraints (such as the LSP being 
uncharged, coupling constants remaining perturbative, no 
Higgs or SUSY particle having been observed at LEP, etc.). 
The extent of these regions is quite limited for the NS mod- 
els, and is largest for the HS models. 

• Second, we determine the masses of the Higgs bosons and 
SUSY particles as a function of (mw2, tan/3). The masses 
of the inos and the sleptons will presumably be measured 
quite accurately, and they will determine the values of 
m 1 / 2 and too; but tan (3 is likely to be poorly determined 
from these masses alone. However, a measurement of m^o 
(or to //o or mjj±) in combination with the mx/2 determi- 
nation from the ino masses will fix a value of tan j3. The 
accuracy of this determination depends upon the accuracy 



with which m^o (and m^o, m#±) can be measured. For 
the A°,H° — > 66 decay modes, for example, this accu- 
racy is fixed by the 66 mass resolution. A resolution of 
±AM;,6 ~ ±10 GeV is probably attainable. For a large 
number, N, of events, to^o can be fixed to a value of order 
AMfcb / V~N, which for N = 20 (our minimal discovery 
criterion) would imply Am^o ^2 — 3 GeV. Such mass 
uncertainty will lead to a rather precise tan (3 determina- 
tion within a given GUT model (except in special cases). 

• Finally, we examine the Higgs branching ratios as a func- 
tion of location in (mi/2,tan/3) parameter space, and de- 
termine the statistical accuracy with which these branching 
ratios can be measured for reasonable assumptions regard- 
ing Higgs tagging and reconstruction efficiencies. 

Ratios of branching ratios are of particular interest since 
certain types of sytematic errors will cancel. Relative Higgs 
branching ratios can be measured by 'tagging' one member 
of the produced pair using a fully reconstructable all-jet decay 
mode, and then looking at the various final states emerging from 
the decay of the other member of the pair. Using the measured 
values of B(h° — > 66) and B(t — > 2jb) and with experimental 
knowledge of efficiencies, we can thus measure 

B(H° ->SUSY)B off (A° -^bb+ti)+B(A" ->SUSY)B cff (if ^bb+tt) 
Soft (if ->bb+tI)B e{f (A ->fcb+tt) 

B(H ^ti)B(A ^bb) + B(A ^ti)B(H°^bb) 
B(H°^bb)B{A°^bb) 

B(H°^h°h )B(A a ^bb) 
B(H°^bb)B{A a ^bb) 

B(A°^Zh a )B(H a ^bb) 
B(H a ^bb)B{A a ^bb) 

B(H + ->SUSY)Bfg~ ->bt)+B(ff- ->SUSY)B(if+ ->tfc) 



B(H+^tb)B(H-^bt) 
B(H + -^r + u)B(H~ ->bt) + B(if~ 



' v)B(H + ->tfe) 



B{H+^tb)B{H- 
B{H + -*h°W + )B(H- ->ht)+B(if" 



->h°W).B(if + ->tfe) 



B(if+->t6)B(if-->6t) 

The SUSY final states can be identified by the presence of miss- 
ing energy opposite the fully reconstructable all jet mode(s) 
used to tag the first member of the Higgs pair. We retain both 
66 and it final states in Eq. ([[]), using an efficiency weighted 
combination denoted by B e g, in order that we may assess the 
importance of SUSY decays both in regions where 66 decays 
of the H°, A are dominant and in regions where it decays are 
important. Note that since to^o ~ m#a we cannot separate the 
H and A decays to the same final state; we can only measure 
the indicated 'average' values of Eqs. ([I]) and (0). 

Two illustrations are provided. In Figs. [I] andg we show con- 
tours in (mi/a, tan/3) parameter space of constant values for 
the ratio of Eqs. ([!]) and (||), respectively, for each of the six 
scenarios defined earlier. The three different curves for each 
value of the ratio indicate the precision with which experiment 
can determine a location in parameter space. These results are 
based on event rates calculated including all relevant branch- 
ing ratios and assuming an 'effective' integrated luminosity of 
L eS = 80 fb _1 at Vs = 1 TcV, where L cS = 80 ftT 1 includes 
an overall tagging, detection, and so forth, efficiency of e = 0.2 



at L = 400 fb _1 (about two years of running). We observe 
that the precision is actually rather good. Since the H°A° and 
H + H~ SUSY ratio contours displayed tend to cross one an- 
other, a measurment of these two ratios will determine a loca- 
tion in (mm, tan/3) parameter space in each GUT scenario. It 
turns out that this determination in one model often disagrees at 
a statistically very significant level with the location determined 
on the basis of the m~± and m<o masses, described earlier, for 
any other model. 



2<B(H°, A°^SUSY)/B eff (H°, A°^bb or tt)> 



e V ^ H°A°. E = I TeV, L„ = 80 fb"' 



c ' 
o 





0.1 






0.6 






2.0 v^ - " 


\ D" 






r 






\a 

















12.5 


(i) 


10 


(2) 


C ™ 

o 

'*-' 5 


(3) 


2.5 


(4) 


40 


(5) 


VL30 


(6) 


B 20 


(7) 


10 




100 150 200 250 300 
HS" 





100 200 300 

m 1/2 (GeV) 



100 200 300 

m 1/2 (GeV) 



Figure 1 : We plot contours along which the ratio of Eq. (jl|) has 
a given constant value, within the constraint/kinematically al- 
lowed (mm) tan j3) parameter space (as indicated by the 'Bdy' 
lines) of the D", D+, NS", NS+, HS", and HS+ models. Re- 
sults are shown for the same three central values for all models. 
For each central value, three lines are drawn. The central line is 
for the central value. The other two lines are contours for which 
the ratio deviates by ±ler statistical error (see Ref. [||]) from the 
central value. Bold lines indicate the boundary beyond which 
fewer than 4 events are found in the final states used to measure 
the numerator of the ratio. 

To more thoroughly illustrate the extent to which the set of 
ratios given in Eqs. (0)-(0) can distinguish between scenarios, 
let us focus on one particular case. Suppose the correct model 
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Figure 2: As in Fig. [IJ, but for the ratio of Eq. (||). 



Table I: We tabulate the values of m 1 / 2 (in GeV) and tan/3 
required in each of our six scenarios in order that rn^a = 
349.7 GeV and m~± = 149.5 GeV. Also given are the cor- 
responding values of m H a and to— . Masses are in GeV. 





D 


D+ 




NS+ 


HS" 


HS+ 


mi/2 
tan/? 


201.7 
7.50 
350.3 
146.7 


174.4 
2.94 
355.8 
127.5 


210.6 
3.24 

353.9 
91.0 


168.2 
2.04 
359.0 
73.9 


203.9 
12.06 
350.1 
222.9 


180.0 
3.83 
353.2 
197.4 



Table II: We give the numbers of events predicted in each 
scenario at the parameter space locations specified in Table | 
available for determining the numerators and denominators of 
Eqs. ([l|)-(0) and Eqs. These event rates are those for 

L c ft = 80 fb _1 at yfs — 1 TeV. They include all branching 
fractions. Our notation is and T>(^ for the event rates 



(#) 



in the numerator and denominator, respectively, of the ratio de- 
fined in Eq. (#). 




D 


D+ 


NS~ 


NS+ 


HS~ 


HS+ 


97.0 


92.3 


88.3 


49.2 


76.1 


124.0 


0.1 


0.7 


3.8 


1.02 


0.0 


0.2 


16.4 


2.7 


46.6 


1.47 


3.8 


2.4 


2.0 


1.3 


9.2 


0.6 


0.4 


1.1 


198 


9.6 


62.1 


2.6 


250 


18.2 


198 


8.9 


58.3 


1.6 


250 


18.0 


225 


189 


138 


135 


189 


262 


58.4 


4.2 


6.5 


1.1 


90.0 


9.5 


13.0 


12.8 


21.9 


9.0 


3.3 


12.3 


317 


415 


445 


465 


320 


348 



is D" with m 1/2 = 201.7 GeV and tan (3 = 7.50. This would 
imply m^o = 349.7 GeV, m~± = 149.5 GeV. The to 1 / 2 and 
tan/3 values required in order to reproduce these same to^o 
and m~± values in the other scenarios are listed in Table Q. Also 
given in this table are the predicted values of mjj° and to~ 
for each scenario. In order to get a first feeling for event num- 
bers and for the errors that might be expected for the ratios of 
interest, we give in Table || the numbers of events, Af and V, 
predicted in each scenario for use in determining the numera- 
tors and denominators of Eqs. and Eqs. assum- 
ing L c g = 80 ftV 1 at Vs = 1 TeV. These numbers include the 
SUSY branching fractions, -B e ff, and so forth. 

In Table ^ we quantify the process of excluding the D + , 
NS~, NS + , HS , and HS + scenarios relative to the input D~ 
scenario. There we give the contribution to A\ 2 (computed rel- 
ative to the assumed-to-be-correct scenario) for each of a se- 
lection of independently measurable ratios. Also given for each 
of the incorrect scenarios is the sum of these contributions. This 
table shows that the D~ scenario can be distinguished from the 
D + , NS~, NS + , and HS + scenarios at an extremely high statis- 
tical level. Further, even though no one of the branching frac- 



tion ratios provides an absolutely clear discrimination between 
the D~ and the HS - scenarios, the accumulated discrimination 
power obtained by considering all the ratios is very substantial. 
In particular, although the ratios of Eq. ([}]), (Q), and (Q) are only 
poorly measured for L c g — 80 ftV 1 , their accumulated A\ 2 
weight can be an important component in determining the like- 
lihood of a given model and thereby ruling out incorrect model 
choices. 

Thus, consistency of all the ratios with one another and with 
the measured to^o, neutralino and chargino masses will gen- 
erally restrict the allowed models to ones that are very closely 
related. The likelihood or probability associated with the best 
fit to all these observables in a model that differs significantly 
from the correct model would be very small. 

An important issue is the extent to which one can be sensi- 
tive to the branching fractions for different types of SUSY de- 
cays of the Higgs bosons, relative to one another and relative to 
the overall SUSY decay branching fraction. Rates in different 
channels depend in a rather detailed fashion upon the SUSY pa- 
rameters and would provide valuable information regarding the 
SUSY scenario. For example, in going from NS to D to HS the 
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Table III: We tabulate A%? (relative to the D _ scenario) for 
the indicated branching fraction ratios as a function of scenario, 
assuming the measured m^o and m~± values are 349.7 GeV 
and 149.5 GeV, respectively. The SUSY channels have been 
resolved into final states involving a fixed number of leptons. 
The error used in calculating each Axf is the approximate la 
error with which the given ratio could be measured for i c ff = 
80 fb _1 at yfs = 1 TeV assuming that the D scenario is the 
correct one. 



Ratio 


D+ 




NS+ 


HS- 


HS+ 


( 




[0£][> 0j]/bb,U 


12878 


1277 


25243 


0.77 


10331 


[li][> 0j]/bb,tt 


13081 


2.41 


5130 


3.6 


4783 


[2£][> 0j]/bb,tt 


4543 


5.12 


92395 


26.6 


116 


h°h°/bb 


109 


1130 


1516 


10.2 


6.2 


H+ 


[oe)[>oj]/tb 


12.2 


36.5 


43.2 


0.04 


0.2 


[ie\[>oj]/tb 


1.5 


0.3 


0.1 


5.6 


0.06 


h°W/tb 


0.8 


0.5 


3.6 


7.3 


0.3 


ru/tb 


43.7 


41.5 


47.7 


13.7 


35.5 


»xf 


30669 


2493 


124379 


68 


15272 



masses of the sneutrinos and sleptons increase relative to those 
for the charginos and neutralinos. The iJ°,^4° — > £ + £~ and 
— ► t^v branching fractions should decline in comparison 
to H°, A — ► xtxl an d H ± — > XiX\> respectively. In small 
sections of the D and NS scenario parameter spaces, the slep- 
tons and sneutrinos are sufficiently light that xt decays almost 



exclusively to £ v followed by £ 



implying that 



xf decays would mainly yield leptons and not jets. 

The difficulty is that several different SUSY channels can 
contribute to any given final state. For example, the £ + £~ +$t 
channel receives contributions from both H°, A — > £ + £~ and 
XiXi decays; and the £^ + JpT channel receives contributions 
from H ± — ► t^v and xtxl- Another example, is the purely 
invisible H or A final state; it can arise from either XiXi or 
vv (with v — > vx\) production. Thus, the physically distinct 
channels, defined by the number of leptons and jets present,[] 
typically have multiple sources. Still, a comparison between the 
rates for the final states so-defined might be quite revealing. For 
instance, if xf —> t^v is not kinematically allowed, the xfx.i 
final states are expected to yield more l£ + 2j and 0£ + 4j events 
than 2£+0j events, whereas £ + £~ events will yield only 2£+0j 
events. Further, the ts must be of the same type in this latter 
case. The effective branching fraction for XiXi ~ * £ + £~+$t 
with both £'$ of the same type is only 1/81. In addition, the ts 
in the latter derive from three-body decays of the xf > an d would 
be much softer on average than ts from £ + £~ . Even if this dif- 
ference is difficult to see directly via distributions, it will lead to 
higher efficiency for picking up the £ + £~ events. 

Based on the above discussion, the following ratios would 



appear to be potentially useful. 

B{H° ->bfe)B(A° -^[0l][0j])+B(A° -^bb)B(H a -> [01] [Oj] ) 
B(ff°-»66)B(A ->SUSY)+B(A°^&6)B(fl" ->SUSY) 

B(H" -^bb)B(A" -^[2e][0]]) + B(A° -^bb)B(H" ->[2l] [0j]) 
B{H°^bb)B(A ^SUSY) + B(A°^bb)B(H ( >^SUSY) 
B(H°^bb)B(A"->[>oe] [0j]) + B(A"^bb)B(H o ^[>ae] [0j]) 
B(H ^bb)B(A ^SUSY) + B(A ^bb)B(H ^SUSY) 

B(H° ->bb)B(A° -*loe\[>lj]) + B(A° -^bb)B(H" ->[0^][>lj]) 
B(H ^bb)B(A ^SVSY) + B(A ^bb)B(H"^ SUSY) 

B(H° -,bb)B(A° '->[!£] [>lj]) + B(A° -^bb)B(H° '^[l^][>lj]) 
B(H ^bb)B(A ^SUSY) + B(A ^bb)B(H a ^SUSY) 

B(H+ ^[U][0j])B(H- ^bt)+B(H- ~,[U][Q 3 ])B(H+ ^ib) 



B(H+ ->SUSY)B(if- ^bt)+B(H~ 
B(H + ^[>U] [0j])B(H- ->bt)+B(H~ 



^SUSY) B{H+^tb) 
,[>U\[0j\)B{H+^ib) 



B(H+^SUSY)B(H-^bt)+B(H- 
B(H + ->[0^][>lj])B(g~ ->bf) + B(g~ 



+SUSY)B(B+— tb) 
>[0l][>lj])B(B/+->t6) 



B(B+^SUSY)B(B-^f)t)+B(if-^SUSY)B(B+^tb) 



(8) 

(9) 
(10) 

(ID 
(12) 
(13) 
(14) 
(15) 



'The totally invisible final state would be [0£] [Oj], and so forth. 



Also of interest are ratios of the different numerator terms to 
one another within the above neutral and charged Higgs boson 
sets. All the ratios that one can form have the potential to pro- 
vide important tests of the Higgs decays to the supersymmet- 
ric particle pair final states. We find that the ratios of rates of 
the various SUSY channels can contribute significantly to our 
ability to discriminate between different GUT scenarios. To il- 
lustrate, we follow the same procedure as in Table III Taking 
myjo = 349.7 GeV and m~± — 149.5 GeV, we assume that 

the correct scenario is D and compute the A% 2 by which the 
prediction for a given ratio in the other scenarios deviates from 
the D~ prediction. Statistics are computed on the basis of the 
expected D~ rates. The resulting A% 2 values are given in Ta- 
ble Since these ratios are not all statistically independent 
of one another, we do not sum their Axf 's to obtain an overall 
discrimination level. However, a rough indication of the level 
at which any given scenario can be ruled out relative to the D~ 
is obtained if we add the largest Axf from the neutral Higgs 
list and the largest from the charged Higgs list. The weakest 
discrimination level following this procedure is Ax 2 ~ 15 in 
the case of the D + scenario. Note that this scenario is highly 
unlikely on the basis of the earlier J^. Axf value listed in Ta- 
ble Ejl In Table Q the weakest discrimination was that for the 
HS~ scenario with ^ i Axf ~ 68. We observe from Table |Tv| 
that the ratio B(H°, A -> [0£][0 j])/B(H°, A -> [2£}[0j]) 
has Axf ~ 928 for the HS~ case, which would certainly rule it 
out. 

The above illustrations demonstrate that the ratios of rates for 
individual SUSY channels correlate strongly with the underly- 
ing physics of the different GUT scenarios (light vs. heavy slep- 
tons in particular) and add a powerful component to our ability 
to determine the correct scenario. 

II. Discussion and Conclusions 

Once the Higgs bosons are detected and their masses deter- 
mined, the relative branching fractions for the decay of a single 
Higgs boson can be measured by 'tagging' (i.e. identifying) one 
member of the H°A° or H + H~ pair in an all-jet mode, and 
then looking at the ratios of the numbers of events in differ- 
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Table IV: We tabulate Ay? (relative to the D~ scenario) for the 
indicated ratios as a function of scenario, assuming the mea- 
sured rriAo and m~± values are 349.7 GeV and 149.5 GeV, re- 
spectively. The SUSY channels have been resolved into final 
states involving a restricted number of leptons and jets. Only 
those ratios with substantial power for discriminating between 
scenarios are tabulated. The error used in calculating each A%f 
is the approximate ler error with which the given ratio could be 
measured for L c g = 80 fb _1 at ^J~s = 1 TeV assuming that the 
D~ scenario is the correct one. 



Ratio 


D+ 


NS~ 


NS+ 


HS- 


HS+ 


(H U ,A U > 


[Of] [Oj] /SUSY 


3.5 


193 


3.4 


1.4 


0.6 


[> 0f][0j]/SUSY 


0.4 


15.3 


6.8 


20.9 


15.8 


[0e][0j]/[2£][0j] 


9.6 


503 


0.1 


928 


105 


[oe\[oj]/[> o£][oj] 


5.8 


41.9 


0.03 


48.4 


24.5 


[Qfl[<y]/[ofl[> u] 


1.4 


1074 


6.4 


3.5 


2.7 


I0l][0j] /[!£][> lj] 


0.3 


3520 


4.3 





1.4 


H+ 


[> li][0j]/SUSY 


1.0 


56.2 


75.2 


3.4 


0.5 


[0£][> ljj/SUSY 


2.1 


21.7 


33.4 


1.3 





[> u][oj]/[oe][> lj] 


5.2 


930 


5738 


4.0 


0.4 



ent event classes on the opposing side. In this way, the relative 
branching ratios of Eqs. (g])-(@), Eqs. Eqs. (|)-(|l|), and 

Eqs. (|l3[)-(|ri|) can be measured with reasonable accuracy when- 
ever parameters are such that the final states in the numerator 
and denominator both have significant event rate. We find that 
the measured Higgs masses and relative branching fractions, in 
combination with direct measurements of the chargino and neu- 
tralino masses, will over-constrain and very strongly limit the 
possible SUSY GUT models. 

The specific SUSY GUT models we considered are moder- 
ately conservative in that they are characterized by universal 
boundary conditions. The strategy for checking the consistency 
of a given GUT hypothesis is straightforward. First, the mea- 
sured A a , neutralino and chargino masses are, in almost all 
cases, already sufficient to determine the mi/a and tan/3 val- 
ues required in the given GUT scenario with good precision. 
The Higgs sector branching fractions can then be predicted and 
become an important testing ground for the consistency of the 
proposed GUT hypothesis as well as for testing the MSSM two- 
doublet Higgs sector structure per se. Typically, a unique model 
among the six rather similar models is singled out by combining 
measurements from the Higgs sector with those from conven- 
tional SUSY pair production. In short, measurements deriving 
from pair production of Higgs particles can have a great impact 
upon our ability to experimentally determine the correct SUSY 
GUT model. 

The above discussion has left aside the fact that for uni- 
versal soft-scalar masses the measured value of the slepton 
mass would determine the relative magnitude of too and TOi/2, 
thereby restricting the possible scenarios (see Table |). How- 

2 In some cases, absolute event rates are so different that they would also pro- 
vide substantial discrimination between different models, despite the possibly 
large systematic errors. 



ever, if the soft-scalar slepton mass is not the same as the soft- 
scalar Higgs field masses at the GUT scale, the branching frac- 
tion ratios would give the best indication of the relative size of 
the soft-scalar Higgs mass as compared to to^/2. 

More information regarding the slepton/sneutrino mass scale 
and additional ability to discrminate between models are both 
realized by subdividing the SUSY decays of the Higgs bosons 
in a way that is sensitive to the relative branching fractions 
for slepton/sneutrino vs. chargino/neutralino decays. Slep- 
ton/sneutrino channels essentially only produce leptons in the 
final state, whereas the jet component is typically larger than the 
leptonic component for chargino/neutralino decays (other than 
the totally invisible x?Xi mode). Thus, we are able to define 
individual SUSY channels, characterized by a certain number 
of leptons and/or jets, which display a strong correlation with 
the slepton/sneutrino decay component. We find that these indi- 
vidual channels have sufficiently large event rates that the ratios 
of the branching fractions for these channels can typically be 
determined with reasonable statistical precision. Excellent dis- 
crimination between models on this basis is found. 

In conclusion, our study shows that not only will detection of 
Higgs pair production in e + e _ or fJ, + /i~ collisions (at planned 
luminosities) be possible for most of the kinematically accessi- 
ble portion of parameter space in a typical GUT model, but also 
the detailed rates for and ratios of different neutral and charged 
Higgs decay final states will very strongly constrain the choice 
of GUT-scale boundary conditions. In estimating experimen- 
tal sensitivity for Higgs pair detection and for measuring Higgs 
masses and branching fractions, we included substantial ineffi- 
ciencies and all relevant branching fractions. Although we be- 
lieve that our estimates are relatively conservative, it will be im- 
portant to re-visit this analysis using a full Monte Carlo detector 
simulation. 
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